function dy=SIR(t,y)
i=y(1),s=y(2),r=y(3);
dy=zer0s(3,1);
dy(1)=y(1)*y(2)-0.1*y(1);
dy(2)=-y(1)*y(2);
dy(3)=0.1*y(1);
